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ABSTRACT 

A variety of observations now provide evidence for vigorous motion in the atmospheres of brown 
dwarfs and directly imaged giant planets. Motivated by these observations, we examine the dynamical 
regime of the circulation in the atmospheres and interiors of these objects. Brown dwarfs rotate rapidly, 
and for plausible wind speeds, the flow at large scales will be rotationally dominated, exhibiting 
geostrophic balance between pressure gradient and Coriolis forces. We present three-dimensional, 
global, anelastic numerical simulations of convection in the interior. Fundamental theory, scaling 
arguments, and our anelastic simulations all demonstrate that, at large scales, the convection aligns 
in the direction parallel to the rotation axis. Convection occurs more efficiently at high latitudes 
than low latitudes, leading to systematic equator-to-pole temperature differences that may reach ~1 
K near the top of the convection zone. The rotation significantly modifies the convective properties. 
The interaction of convection with the overlying, stably stratified atmosphere will generate a wealth 
of atmospheric waves, and we argue that, just as in the stratospheres of planets in the solar system, 
the interaction of these waves with the mean fiow will lead to a significant atmospheric circulation at 
regional to global scales. At scales exceeding thousands of km, this should consist of geostrophically 
balanced, stratified turbulence (possibly organizing into coherent structures such as vortices and jets) 
and an accompanying overturning circulation. We present a semi-quantitative, analytic theory of this 
circulation as a function of the wave-driving efficiency. For plausible efficiencies, this theory predicts 
characteristic horizontal temperature variations of several to ^50 K, horizontal wind speeds of ^10- 
300 ms~^, and vertical velocities that advect air over a scale height in ~10^-10^ s. The vertical mixing 
associated with this large-scale circulation may help to explain the chemical quenching of CO and 
NII3 observed on some brown dwarfs. Moreover, the implied large-scale organization of temperature 
perturbations and vertical velocities suggests that, near the L/T transition, patchy clouds can form 
near the photosphere, helping to explain recent observations of brown-dwarf variability in the near- 
infrared. 

Subject headings: stars: low-mass, brown dwarfs; convection; waves; turbulence 
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1. INTRODUCTION 

Since the discovery of brown dwarfs beginning in the 
mid-1990s, our understanding of the atmospheric struc- 
ture of these objects has grown ever more sophisticated. 
Approximately 1000 brown dwarfs, and a handful of di- 
rectly imaged planets, have now been discovered. Obser- 
vational acquisition of infrared (IR) spectra for many of 
these objects have allo wed the de finition o f the L, T, and 
Y spe ctral classes (e.g., Kirkpatrick 2005 j Gushing et al. 
2011 1. The theory for these objects now encompasses a 
broad understanding of their evolution, radii, luminosity, 
molecular composition, spectra, and colors, and includes 
prescriptions for condensate formation and rainout, sur- 
face patchiness, and disequilibrium chemistry. Notably, 
however, these theoretical advances have relied heavily 
on one-dimensional (ID) models for the atmospheric ra- 
diative tr ansfer and interior ev olution (for an early re- 
Burrows et al.||2001|). By comparison, little 



view, see 



effort has been made to understand the global, three- 
dimensional atmospheric dynamics of these substellar 
bodies. 

Yet there is increasing evidence that brown dwarfs ex- 
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hibit vigorous atmospheric circulations. This evidence 
falls into three main classes. First, L dwarfs, particularly 
of later spectral type, show a reddening of near-infrared 
(e.g., J — K) colors that indicate the presence of silicate 
clouds in the visible atmosp heres (e.g ., Kirkpatrick et al.| 
[1999; Kirkpatrick 2005; Cha brier et al. 2000; Tsuji_200: 



tlushing et al. 2006; Knapp et al.||2004p. Since 'clon 



particles would gravitationaily settle in the absence of 
dynamics, such clouds imply the presence of atmospheric 
vertical mixing necessary to keep the particles suspended. 
In the cooler T dwarfs, the condensation occurs progres- 
sively deeper and, for objects with sufficiently low ef- 
fective temperature, eventually no longer influences the 
infrared spectrum. However, the L/T transition itself 
remains poorly understood; it occurs over a surprisingly 
small range of effective temperature and accompanies a 
J-band brightening tha t are not easily captured by stan- 
CChabrier et al.|2000i lAllard et ai.|2001 
06b||Saumon &: Marley 2008 ). Hypothe- 



dard ID models 



Burrows et al.|2006b||Saumon &: MaiTey 2008). Hypothe 
ses that have been put forward to resolve this discrep- 
ancy are that, across the transition, the clo ud sedimen- 
tation efficiency changes ( Knapp et al.||2004 ) or that the 
clouds become patchy, allowing contributions from both 
cloudy and cloud-fr ee regions to affect the disk-integrated 
emergent spectrum ( Ackerman fc Marley 2001 Burgasser 
et al.||2002 [Marley et al.||2010[ ). in both CcLSGS, cL role tor 



atmospheric dynamics m modulating the clouds is impli 
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2 



Showman and Kaspi 



A second line of evidence for atmospheric circulation 
comes from chemical disequilibrium of CO, CH4, and 
NH3 inferred for many cool brown dwarfs. Late T 
dwarfs have sufhciently cool atmospheres that the pre- 
ferred chemical-equilibrium forms of carbon and nitrogen 
near the photosphere are CH4 and NH3, respectively; in 
contrast, CO and N2 dominate under the high-pressure 
and temperature regions at depth. Fitting of infrared 
(IR) spectra to radiative transfer models shows that, near 
the photosphere, many T dwarfs exhibit an overabun- 
dance of CO and an underabundance of NH3 relative to 
chemical equilibrium. This can be attributed to verti- 
cal transport of CO-rich and NHa-poor air from depth 
and the subsequent chemical quenching of these disequi- 
librium mixing ratios due to the long chemical intercon- 
version timescales in the low-pressure, low-temperature 
regions near the photosp here. This sto ry was first worked 
out for CO on Jupiter d P rinn fc Ba r"shay^^l977, Bezard, 



et al.||2002[ |Visscher fc Moses 2U11) and tEen for both 
UUTnd JNH.s on Gl m h (Fcglc y fc Lodders|[l996ljNdl 
eT^pQTl iGrifhth fc Y cUc 1990| [Saumon et al.p(300[ ). 
yubsequently, chemical disequilibrium and vertical mix 



ing have bee n inferred in the atmosphere s of a wide range 
of T dwarfs dSaumon et al. 2006, 2007: Hube ny fc Bur- 
][2007| [Leggett et al.|,2007b,a, ,2008, ,2010t |l>tephe~ 



rows 
et al 



2009[ ). 

i'hird, ^recent near-IR photometric observations 
demonstrate that several brown dwarfs near the L/T 
transition exhibit large-amplitude variability over rota- 
tional timescales, probably due to cloudy and relatively 
cloud-free patches rotating in and out of view. The pos- 
sibility of weather on brown dwar fs has long mo t ivated 
searches for variability. Recently, Artigau et al. (20091 
observed the T2.5 dwarf SIMP0136 in J and bandlind 
found peak-to-peak modulations of ~5% (^5 mmag) 
throu g hout the inferred 2.4-hr rotation period. [Radigan] 
et al I ( |2012[ ) observed the T1.5 dwarf 2M2139 in J, H, 
and Kg and found peak-to-peak variations of up to ~25% 
with an inferred rotation period of either 7.7 or 15.4 hr. 
The relative amplitudes of the variability at J, H, and Ks 
place strong constraints on the cloud and thermal struc- 
ture associated with the variability. These authors con- 
sidered models where the variability resulted from lateral 
variations in effective temperature alone (with no vari- 
ations in the cloud properties), lateral variations in the 
cloud properties alone (with no variation in effective tem- 
perature) , and lateral variations in both temperature and 
cloud properties. The observations rule out models with 
a uniform cloud deck and instead strongly favor models 
with significant lateral variations in both cloud opacity 
and effective temperature; the relatively cloud-free re- 
gions exhibit effective temperatures ~ 100-400 K greater 
than the cloudier regions. This suggests a picture with 
spatially distinct regions of lower and higher condensate 
opacity, where radiation escapes to space from lower- 
pressure, cooler levels in the high-opacity regions and 
deeper, warme r levels in the low-op acit y regions . The 



observations of Artigau et al. ( 2009 ) and Radigan et al. 
( [2012) both show that the lightcurves vary significantly 
over mtervals of several Earth days, indicating that the 
shape, orientation, or relative positions of the low- and 
high- condensate opacity regions evolve over timescales of 
days. 

In addition to these observations of field brown dwarfs. 



growing numbers of young, hot extrasolar giant plan- 
ets (EGPs) are being imaged and characterized. Promi- 
nent discoveries include planetary-mass companions to 
/3 Pic, 2M1207, and HR 8799. Mufii-band photome- 
try already indicates that 2M1207b and several of the 
HR 8799 planets exhibit clouds and probably disequilib- 



(Hinz et al.||2010 |Bowler et al.||2010 |Currie et al.||2011 


Caficher et al.| 20111 Skemer et al.||2011 
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et al. 2011| Barman et al. 2011a|b |Mar 


ey et al. 2012p. 



planets radiate IR fluxes > 10^ Wm~^, orders of magni- 
tude greater than the flux received by their primary star. 
Stellar irradiation is therefore negligible to their dynam- 
ics. From a meteorology perspective, this population of 
bodies will therefore resemble low-mass, low-gravity ver- 
sions of free-floating brown dwarfs. With next-generation 
telescope facilities, including the Gemini Planet Imager 
and SPHERE, significant numbers of new planets will 
be discovered, greatly opening our ability to probe plan- 
etary meteorology at the outer edge of stellar systems. 

These existing and upcoming observations provide 
strong motivation for investigating the global atmo- 
spheric dynamics of brown dwarfs and directly imaged 
planets. As yet, however, no investigations of the global 
atmospheric circulation of brown dwarfs have been per- 
formed. The only study of brown dwarf atmospheric dy - 
namics published to date is that of [Freytag et al. (2010), 



who performed two-dimensional, non-rotating convection 
simulations in a box ^400 km wide by ~150km tall. 
Their study provides valuable insights into the role of 
convectively generated small-scale gravity waves in caus- 
ing vertical mixing. Nevertheless, dynamics on scales of 
tens to hundreds of km differs substantially from that on 
global scales of 10'*-10^ km, and thus, for understand- 
ing the global-scale circulation — including the implica- 
tions for variability — it is essential to consider global- 
scale models. 

Here, we aim to fill this gap by presenting the first 
global-scale models of brown dwarf atmospheric dynam- 
ics. Rotation periods of L and T dwarfs inferred from 
spectral line broadeni ng range fr om ^2 to 12 hours ( Za- 



patero Osorio et al. 2006; Reiners fc Basri 20081, m 



fine with the periods oT SIM P0136 and 2iVL213glnferred 



from lightcurve modulation ( Artigau et al. 2009 Radi- 
We will show that, at these periods, 



gan et al. 



2012) 



rotation dominates the global-scale dynamics and will 
constitute an overriding factor in controlling the three- 
dimensional wind and temperature structure. We first 
present theoretical arguments to highlight the fundamen- 
tal dynamical regime in which brown dwarfs lie, to show 
how rotation organizes the large-scale wind, establishes 
systematic temperature differences, and shapes the con- 
vective velocities (Section [2]). We next present global, 
three-dimensional numerical simulations of the interior 
convection of brown dwarfs that confirm our theoretical 
arguments and provide insights into the detailed, time- 
evolving global wind and temperature patterns and dy- 
namical timescales (Sectionjsl). We then consider the dy- 
namics of the stably stratined atmosphere that overlies 
the convective interior, demonstrating how large-scale 
vortices and/or zonaj^jets are likely to emerge from in- 

^ Zonal refers to the east-west (longitudinal) direction whereas 
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teractions with the interior (Section |4]). Next, we con- 
sider observational imphcations, since IR radiation to 
space typicaUy occurs from within this stratified layer 
and hence infrared spectra and lightcurves are strongly 
shaped by its dynamics (Section pi). The final section 
concludes (Section [6]). 

We emphasize that our goal is to provide a theoretical 
foundation for understanding the atmospheric dynamics 
of rapidly rotating, ultracool dwarfs and young EGPs, 
broadly defined. As such, we emphasize dynamical con- 
siderations, and intentionally simplify our models by ex- 
cluding clouds, chemistry, and detailed representation of 
radiative transfer. This provides a clean environment in 
which to identify key dynamical processes and construct 
a theoretical foundation for more realistic studies that 
will surely follow. 

2. BACKGROUND THEORY: APPLICATION TO 
BROWN DWARFS 

Here, we review basic concepts in atmospheric fiuid 
dynamics and apply them to brown dwarfs to understand 
the large-scale structure of the fiow. 

2.1. Importance of rotation 

Brown dwarfs rotate rapidly, and this will exert a ma- 
jor influence on their atmospheric dynamics. To demon- 
strate, consider the momentum equation for a rotating 
fluid, which is given in the rotating reference frame of 
the brown dwarf by 



^ +2rj X v = --Vp- V*, 
at p 



(1) 



where dv/dt = dv/dt + v ■ Vv is the material deriva- 
tive, V is the 3D velocity vector, 17 is the planetary ro- 
tation vector, <& is a force potential which includes both 
the gravitational and centrifugal accelerations, V is the 
three-dimensional gradient operator, t is time, and p and 
p are the p ressure and density fields respectively (Ped- 
losky 1987). For the purpose of this discussion we will 
assume the dynamics are inviscid. We also for the present 
neglect the Lorentz force, as appropriate for the atmo- 
spheres and molecular envelopes of cool brown dwarfs; 
nevertheless, the Lorentz force will be crucial in th e dee p 
interior, and return to a discussion of it in Section [XS] 

The nature of the flow depends on the Rossby num- 
ber, given by the ratio of the advective and Coriolis 
accelerations, Ro — U /^L, where U is a. characteris- 
tic wind speed, L is a characteristic length scale of the 
flow, and Q is the rotation rate {2tt over the rotation pe- 
riod). If Ro <^ 1, the flow is rotationally dominated; if 
Ro ^ 1, rotation is important but not dominant, while 
if Ro ^ 1, rotation plays a minor role (see, e.g., Vallis 
2006[ p. 84). Because of their fast rotation rates, Jupiter, 
Saturn, Uranus, and Neptune, as well as the extratropi- 
cal atmospheres of Earth and Mars, all exhibit Ro <^ 1 
at large scales; the only solar-system atmospheres where 
rotation does not dominate are Titan and Venus. 

Estimating Ro for brown dwarfs requires knowledge 
of wind speeds and flow length scales, which are un- 

meridional refers to the north-south (latitudinal) direction. Zonal 
and meridional winds are winds in the eastward and northward 
directions, respectively. Zonal jets refer to atmospheric jet streams 
oriented in the east-west direction. 
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Fig. 1. — Rossby numbers expected on brown dwarfs as a function 
of rotation period and characteristic wind speed, assuming dynam- 
ical features are global in scale (length equals one Jupiter radius). 
Rossby numbers range from ~10~* to 0.1, indicating that the 
regional- and global-scale dynamics in brown-dwarf atmospheres 
will be rotationally dominated over a wide range of parameters. 



knownj^ Nevertheless, Artigau et al. (2009) show that 
if the evolution of their light curve shapes over intervals 
of days is interpreted as differential zonal advection of 
quasi-static features by a latitude-dependent zonal wind, 
the implied differential rotation is ^1%, which for the 
2.4-hour rotation period of SIMP0136 implies a zonal 
wind speed of ~300-500ms~^ depending on the latitude 
of the features. A similar analysis by jRadigan et al.| 
(2012 1 suggests a possible zonal wind speed of ~45ms~^ 



for 2M2139, although they caution that this estimate re- 
lies on rather tentative assumptions. Interestingly, these 
values bracket the range of wind speeds measured for 
the giant planets in the solar system, which range from 
typical speeds of ^30m s~^ o n Jupiter to ~300ms~^ on 
Neptune (e.g., Ingersoll||l990[ )p] Later, we show that the 
large-scale winds in the convective inte rior of a brown 
dwarf are likely to be weak (Section |3.3[ ), but that winds 
potentially exceeding ~10^ms^^ could develop in the 
stratified atmosphere (Section B. 

Regarding length scale, the fact that SIMP0136 and 
2M2139 exhibit large-amplitude variability hints that at- 
mospheric features could be near-global in size (partic- 

* Here, we seek to understand the global-scale flow and the 
appropriate values are therefore not the convective velocities and 
length scales but the wind speeds and length scales associated with 
any organized jets (i.e. zonal flows) and vortices that may exist. 

^ The maximum observed speeds, expressed as a difference be- 
tween the peak eastward and peak westward zonal winds, are sev- 
eral times these typical values, reaching '~200ms~^ on Jupiter and 
~600ms~^ on Neptune. 
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ularly for 2M2139 where variabihty reaches 25%). This 
would miply L ^ Rj ^ 7 x 10^ m, where Rj is Jupiter's 
radius. On the other hand, only a small fraction of 
brown dwarfs exhibit such large variability and it is pos- 
sible that length scales are typically smaller; for example, 
L ^ 10^ m on Jupiter, Saturn, Uranus, and Neptune. 

Adopting length scales L ~ 10^-10^ m, wind speeds 
U ^ 10-1000 ms~^, and rotation rates of 2 to 10 hours 
yields Rossby numbers ranging from 0.0001 to 0.4. Fig- 
ure [T] shows the Rossby number as a function of wind 
speed and rotation period for the case of global-scale 
flows. The values are much less than one everywhere 
except for the largest wind speeds and slowest rotation 
periods considered. This implies that, in general, the 
large-scale circulation on brown dwarfs will be rotation- 
ally dominated. The Ro <C 1 condition on brown dwarfs 
implies that the flow is geostrophically balanced, that 
is, the primary balance in the momentum e quation is 
betwe en Coriolis and pressure-gradient forces (Pedlosky 
T987|. 



2.2. Organization of flow 

Significant insight into the flow structure can be ob- 
tained from the vorticity balance. Taking t he curl of 
Tl gives a vorticity equation of the form (e.g., Pedlosky 



'87 



In 



•Vv-(20 



Vp X Vp 



(2) 



where w = V x v is the relative vorticity. The term 
on the right side, called the baroclinic term, is nonzero 
when density varies on constant-pressure surfaces. Note 
that, since w scales as C//L, the ratio w/fi ~ Ro, and 
the time-derivative term is generally order Ro smaller in 
magnitude than the second and third terms on the left 
side. Taking Ro <^ 1, appropriate to the flow on a brown 
dwarf, yields a leading-order vorticity balance given by 



Vp X Vp 



(3) 



One might expect that convection homogenizes the en- 
tropy within the convection zone, in which case density 
does not vary on isobars and Vp x Vp = 0. This is 
called a barotropic flow. In this case, Eq. ([3 1 simply be- 
comes the compressible- fluid generalization of the Taylor- 
Proudman theorem, which, expressed in a cylindrical co- 
ordinate system centered on the rotation axis, is 



du 
dz 



dz 



= 



(4) 

• = 0. (5) 

where z is the direction parallel to the rotation axis, 
u is the azimuthal (zonal) velocity, v± the velocity to- 
ward/away from the rotation axis, vj^ — (u, v±) is the ve- 
locity in the plane perpendicular to the rotation axis, and 
V_L is the gradient operator in the plane perpendicular to 
the rotation axis. The theorem states that if the flow has 
a small Rossby number and is invicsid and barotropic the 
fluid motion will be completely two-dimensional, there- 
fore there will be no variation in the fluid velocity along 
the direction of the axis of rotation. The fluid then moves 
in columns aligned with the rotation axis. Note that no 



constant-density assumption was made; Equations Q- 
^ hold in a barotropic, geostrophic, low-viscosity fluid 
even if the density varies by orders of magnitude across 
the system. 

Within a spherical planet or brown dwarf, such 
columns can easily move in the zonal (east-west) direc- 
tion. However, the columns cannot easily move toward 
or away from the rotation axis, because this changes 
the length of the columns and the local density within 
them, both of which induce non-zero Vj_ • vj^ that vio- 
late the theorem. For such a barotropic fluid, the pre- 
dominant planetary-scale circulation therefore consists of 
zonal (east-west) wind whose speed varies minimally in 
the direction along the rotation axis; by comparison, the 
north-south flow is weak. Of course, the theorem is only 
valid to order Ro, and so motions toward/away from the 
rotation axis — as well as motions parallel to the fluid 
columns — can occur but only with amplitudes ^Ro less 
than that of the primary zonal flow|^ 

In reality, turbulent convection results in horizontal 
entropy gradients and therefore the fluid is not in a 
barotropic state leading to a non-vanishing term on the 
right side of Eq. ([3|. As a result, shear can develop 
along the z direction. Considering the zonal component 
of Eq. ^ yields 



2n^ = -^l^.x 

oz 



(6) 



where A is the unit vector in the longitudinal direction. If 
the flow exhibits minimal variation in longitude, then it 
can be shown that |VpxVp| = |Vp|(9p/9y)p. Since |Vp| 
is overwhelmingly dominated by the hy drostatic compo- 
nent, we have to good approximation ( [Showman et al^ 
2010t 

dz pr \dcj)J p 

where (j) is latitude and r is radial distance from the cen- 
ter of the planet. Thus, variations in the geostrophic 
wind along z must be accompanied by variations in den- 
sity on isobars. This relation, well-known in atmospheric 
dynamics, is called the thermal-wind equation. 

By itself, however, the preceding theory gives little in- 
sight into the spatial organization — columnar or not — of 
the internal entropy perturbations and any thermal- wind 
shear that accompanies them. An alternative point of 
view that sheds light on this issue is to consider the an- 
gular momentum budget. The angular momentum per 
unit mass about the rotation axis is given by 



M = A'U 



'J-n 



M, 



= Qr^ cos^ (j) + ur cos cf). 



(8) 



where the flrst and second terms represent the contri- 
butions due to the planetary rotation and winds in that 
rotating frame, respectively. Writing the zonal momen- 
tum equati on in terms of angular momentum yields (e.g., 
Peixoto fc Oort„1992( Chapter 11) 

p"-^--'!- (9) 



dt 



^ Of course, at very small scales, the Rossby number exceeds 
unity and the convection at these sma ll scales will not organize into 
columns (e.g., |Glatzmaier et al.|2009| . The columnar organization 
applies only at scales suHiciently large that Ro <C 1- 
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It is useful to decompose the pressure and density into 
contributions from a static, wind-free reference state and 
the deviations from that state due to dynamics. When 
wind speeds are less than the speed of sound, these dy- 
namical density and pressure perturbations are small, 
leading to a continuity equation V • (pv) = 0, where p is 
the reference density profile[^ Motivated by the fact that 
the convective eddies drive a mean flow, we represent 
the dynamical variables as the sum of their zonal means 
(denoted by overbars) and the deviations therefrom (de- 
noted by primes), such that M = M+M' , v = v-|-v', etc. 
Here, we refer to these overbarred quantities as the mean 
flow and the primed quantities as the eddies. Substitut- 
ing these expressions into Eq. ([9| and zonal ave raging 
leads to th e zonal-mean momentum equation (cf Kaspi 
eral][2009| 



But how important are the eddy torques on the right- 
hand side of Eq. (12)? At small Rossby number, the 



dM 



+ V • (pv A/) + V • (pv'M') = 0. 



(10) 



This equation states that temporal changes to the zonal- 
mean angular momentum at any given location (first 
term) can only result from advection of the zonal-mean 
angular momentum by the mean flow (second term) or 
changes in the zonal-mean angular momentum due to 
torques caused by eddy motions (third term). In a sta- 
tistical steady state, M equilibrates and the equation 
simply becomes 

V • (pvM) + V • (pVIF) = 0. (11) 

Since the ratio of M„ to Mfi is essentially the Rossby 
number, it follows that for the Ro <C 1 regime expected 
on a brown dwarf, M ^ Mq. Thus, for a rapidly rotating 
brown dwarf, surfaces of constant angular momentum are 
nearly parallel to the axis of rotation. Using this result, 
along with the continuity equation, Eq. (Ill becomes 



V- VAfn = - ^V • (pv'M'). 

P 



(12) 



This result has major implications for the circulation 
on rapidly rotating giant planets and brown dwarfs. The 
equation implies that the mean flow, v, can only cross an- 
gular momentum surfaces in the presence of eddy corre- 
lations between v' and M' (such eddy correlations cause 
a torque that changes the zonal-mean angular momen- 
tum following the mean flow, as necessary for the mean 
flow to cross angular momentum surfaces). When such 
eddy effects are small, or if the flow is axisymmetric with 
no variation in longitude (for which v' and M' are zero 
by definition), then v- VMn = 0. In such a situation, the 
mean flow must to leading orde r be parallel to surfaces 
of con stant angular momentum (lBusse|1976 Kaspi et al. 
20091 and there can be no flow crossing these surfaces, 
i'his constraint places no limitation on the zonal-mean 
zonal flow, u, but requires the meridional circulation v 
to be small. It is important to emphasize that this con- 
straint differs from the Taylor-Proudman theorem, since 
it does not require the flow to be barotropic, nor does it 
state that u is independent of z. 

The reference density will generally be a function of radius. 
Note that this approximate continuity equation (essentially the 
anelastic approximation) filters sound waves from the system, 
which is a reasonable approximation as long as wind speeds are 
less than the speed of sound. 



planetary rotation contains so much angular momentum 
that even in the presence of vigorous convection, eddy 
torques are unable to drive a rapid mean-meridional cir- 
culation; therefore, we still have v • VMq « at leading 
order. To show this, we can estimate the timescale of the 
meridional circulation and compare it to the characteris- 
tic timescale for convection. The characteristic timescale 
for convection to traverse the interior is Tconv ~ -D/w, 
where Z? is a thickness of the layer in question (e.g., the 
planetary radius) and w is the characteristic convective 
speed. We can estimate the timescale for the mean flow 
to cross angular momentum contours as follows. The 
eddy correlation v'M' is just v'u'r_L, where rj^ is the dis- 
tance from the rotation axis. Under the assumption that 
the eddy velocities scale with the convective velocities 
w, we can write v'M' ^ Cr ^uP' where C is a correla- 
tion coefficient equal to one when u' and v' are perfectly 
correlated and equal to zero when v! and v' e xhibit no 
correlation. To order of magnitude, Eq. (12 1 then be- 
comes 

v^^r^ ~ Cw^ (13) 

and the meridional velocity therefore has a characteristic 
magnitude 

Defining a timescale for the meridional circulation, 
Tmcrid = T^jv^-, Implies that 



^merid ^ "^ci 



which can be expressed as 



'^mcrid 



CRoc 



(15) 



(16) 



where Rocom — w/rj_D, is a convective Rossby number 
giving the ratio of the convective velocities to the typ- 
ical rotational velocity of the planet in inertial space. 
For typical brown dwarfs, where rotational velocities are 
tens of kms~^, we expect Roconv ^ 1; given the expected 
convective velocities (see we expect Roconv ~ 10~^- 
10~^. Therefore, the timescale for the meridional cir- 
culation is several orders of magnitude longer than the 
characteristic convection timescale. 

This also means that the convective heat transport will 
be more efficient along (rather than across) surfaces of 
constant angular momentum. Fig. [2] shows the onset of 
convection for an experiment driven by a constant heat 
flux at the bottom boundary. Two models are shown, a 
rapidly rotating case on the left and a slowly rotating case 
on the right. While for the larger Rossby number case 
the dominant driving force for the turbulent plumes is the 
buoyancy and therefore the plumes are driven away from 
the center of gravity, for the small Rossby number case 
the convection becomes aligned along the direction of the 
axis of rotating demonstrating the angular momentum 
constraint (v • VMo — 0). Note that for small Rossby 
number experiment only close to the boundaries does 
the Rossby number approach one and therefore there the 
convective cells can close. Thus rotation strongly mod- 
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ulates the heat transport from the interior of the brown 
dwarf at large scales. 



3. THREE-DIMENSIONAL CIRCULATION MODEL 
OF CONVECTION-ZONE DYNAMICS 

3.1. Model 

We solve the fluid equations for a self-gravitating, con- 
vecting, three-dimensional, rota ting bro wn dwarf. We 
adopt the anelastic system (e.g., Ogura fc PhiUips 1962; 
Gilman fc Glatzmaier||1981[ |Ingersoll fc' Poiiard„1982| ), 
which assumes that dynamics introduces only small per- 
turbations of the density, entropy, and pressure from a 
specified reference state, which we here take to be isen- 
tropic. Dynamical density perturbations then enter the 
momentum equations in the buoyancy term but do not 
appear in the continuity equation; this has the effect of 
filtering acoustic waves from the system. The anelastic 
system is appropriate for the fluid interior of a brown 
dwarf, where dynamical perturbations of entropy, den- 
sity, and pressure due to convection are expected to be 
modest and convection should lead to a nearly constant 
entropy throughout. Although this study represents its 
first application to brown dwarfs, the anelastic system 
has previously been used with great success for under- 



standing convection in Jupiter and Saturn (Kaspi et al 
2009] [Jones fc Kuzanyan||2009t |Glatzmaier et al. 2009 
Showm an et al. 2011) and stellar interiors ( jMiesch fc, 

I'oomre 2009, and references therein). 

Our particular implementation is that of |Kaspi et al.| 
(|2009). The momentum, continuity, and energy equa- 
tions, respectively, are given by 

c) 1 ' 1 

-^ + (2fJ+a;)xv = --Vp'-^V^--Vv^+iyV'^v (17) 

ot p p 2 



V • (pv) = 

ds' 1„ „ 1^ 
ot p p 



(P^V,s') = I 



(18) 
(19) 



where Q is thermodynamic heating/cooling per mass, v 
is the kinematic viscosity, k is the thermal diffusivity, and 
other quantities are as defined previously. Here, both v 
and K are taken as constants and are intended to parame- 
terize small-scale eddy mixing. The quantities p(r), p{r), 
and T(r) are the radially varying reference profiles of 
density, pressure, and temperature, respectively; p' and 
p' are the deviations of the density and pressure from 
their local reference values, such that the total pressure 
and density are p = p + p' and p = p + p' . Likewise, s' is 
the deviation of entropy from its reference state value. 

The system is closed with an equation of state, which 
enters through the reference profiles and through the 
relationship between the density, entropy, and pressure 
perturbations in the anelastic system 



P_ 
P 



1 /dp 
P \ds 



dp 
dp 



P 



+ I3p', (20) 



where as and j3 are the radially varying isobaric entropy 
expansion coefficient and isentropic compressibility, re- 
spectively, along the model's radially varying reference 



profile. Here, we adopt the SCVH e quation of state 
(EOS ) for hydrogen-helium mixtures (Saumon et al. 
19951). Given a specified internal entropy and the as- 



sumption that the reference state is in hydrostatic bal- 
anc^ this EOS allows us to calculate t he radially vary- 
ing reference profiles p, p , and T (e.g., Guillot fc Morel 
1995 Guillot et al.||2004 ) as well as the radial profiles of 
Qfs and p along the reference adiabat. 

Many studies of convection in rotating spherical shells 
force the system by passing a heat flux through imperme- 
able upper and lower boundaries, either with a constant- 
temp erature or constan t heat flux boundary condition 



Christensen' l2001[ |2002| lAurnou fc 6lson||2001 

I, and 



Heimpel et al. 2005, and many others) . However, this is 



unrealistic in the context of a substellar object. At high 
Rayleigh numbers, passing a heat flux through the model 
boundaries will lead to thin hot and cold boundary lay- 
ers at the bottom and top boundaries, respectively, which 
detach and form hot and cold convective plumes that in 
some cases can dominate the dynamics. Because real 
brown dwarfs are fluid throughout, the bottom bound- 
ary layer, in particular, is unrealistic. Instead, we force 
the system by imposing a vertically distributed source 
of internal heating and cooling throughout the bottom 
and top portions of the domain, respectively, thus al- 
lowing outward convective transport of heat without the 
development of artiflcial boundary layers (for more de- 
tail see ,Kaspi et al.,2009) p] The top and bottom thermal 
boundary conditions correspond to zero heat flux. The 
top and bottom mechanical boundary conditions are im- 
permeable in radial velocity and free-slip in horizontal 
velocity. 

We solve the equations in spherical geomet ry using 
the state-of -the-art circula t ion m odel MITgcm (Adcroft 
et al. 2004), which Kaspi ( 2008| adapted for anelastic 
simulations of the deep convective envelopes of giant 
planets. The equations are solved using a finit e- volume 
discretization on a staggered Arakawa C grid ( Arakawa 
fc Lamb|1977 ) in longitude and latitude. Our typical res- 
olution is 1" in longitude and latitude with 120 vertical 
levels spaced to give enhanced resolution near the top of 
the domain where t he pressure and de nsity scale heights 
are the smallest (see Kaspi et al.|2009 ). Most models ex- 
tend the full 360° in longitude and in latitude from 80° S 
to 80°N. For some parameter variations, we performed 
simulations in sectors 90° of longitude wide (using a pe- 
riodic boundary condition between the eastern and west- 
ern boundaries) with a resolution of 2° of longitude and 
latitude and 120 vertical levels. For all models, the outer 

* We emphasize that the dynamical model itself is non- 
hydrostatic; hydrostatic balance is used only in defining the ref- 
erence state. 

^ Real brown dwarfs of course do not have substantial internal 
heat sources (the burning of deuterium not playing a role except 
in the most massive objects) but rather decrease in entropy very 
gradually over multi-Byr timescales. Thus, ideally, one would like 
to set up the problem with a heat sink near the top (i.e., cooling) 
and no heat source near the bottom, thereby allowing the internal 
entropy to decline with time in a brown-dwarf like fashion. The dif- 
ficulty is that, due to computational limita tions, achieving steady 
state requires the system to be overforced ( [Showman et al |2011| , 
and without a source of energy near the bottom the interior en- 
tropy would decline unrealistically rapidly. Adding a heat source 
near the bottom, as we have done, allows the global-mean inte- 
rior entropy to be essentially constant over dynamical timescales, 
consistent with expectations for brown dwarfs. 
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and inner boundaries are spherical surfaces with radii of 
IRj and 0.5i?j respectively. This choice of inner bound- 
ary is sufficiently deep to minimize any artificial effect of 
the lower boundary on the surface dynamics. We gen- 
erally use a Jupiter-like interior reference profile with a 
pressure at the outer surface of 1 bar and 20 Mbar at the 
bottom boundary. The interior reference density varies 
by a factor of over 10"^ from t he 1-bar level to the deep 
interior (see Kaspi et al.|2009 ). All simulations are spun 
up from rest using an initial thermal profile correspond- 
ing to the reference profile, and are integrated until a 
statistical steady state is achieved. 

3.2. Results: Convective and thermal structure 

We perform simulations using rotation periods ranging 
from 3 to 200 hours (spanning the typical range observed 
for brown dwarfs) as well as an additional sequence of 
parameter variations adopting rotation periods as long 
as 2000 hours to illustrate the effect of rotation on the 
dynamics. 

Before presenting models for fully equilibrated brown 
dwarfs, we first demonstrate with a pedagogical example 
the crucial importance that rotation plays in the brown- 
dwarf parameter regime. Figure [2] depicts the tempera- 
ture structure during the spin-up phase for two models 
that are identical except for the rotation period, which 
is 2000 hours in the model on the right and 10 hours for 
the model on the left. The models in Figure [2] are not in- 
tended to be realistic brown dwarf models (for example, 
they are forced by a heat flux from the bottom boundary, 
which is not realistic in the context of a brown dwarfs) 
but are instead simply an illustration of the importance 
of rotation in the brown dwarf parameter regime. Nev- 
ertheless, the models do have a realistic Jovian interior 
structure, with density increasing by a factor of ^10^ 
from the interior to the exterior. 

In the slowly rotating model (right panels of Figure^, 
the Coriolis forces are sufficiently weak that the Rossby 
number is >1, so that rotation plays a negligible role in 
the dynamics. Convective plumes rise from the lower 
boundary and ascend quasi-radially toward the outer 
boundary. The plumes are equally able to radially tra- 
verse the domain whether they emanate from the polar 
or equatorial regions, and to zeroth order, the convection 
appears to be isotropic. 

In contrast, the rapidly rotating model exhibits 
Ro <^ 1 and is thus rotationally dominated (Fig.^ left 
panels). As predicted by the theory in Section^ the 
convection develops a columnar structure. Plumes that 
emerge in the polar regions can ascend and descend ra- 
dially while remaining at nearly constant distance from 
the rotation axis; they therefore easily traverse the do- 
main. But plumes forming at lower latitudes cannot eas- 
ily cross the domain because doing so requires them to 
change distance significantly from the rotation axis. In 
a Ro <^ 1 flow, angular momentum is not homogenized 
and lines of constant angular momentum are nearly par- 
allel to the rotation axis. As a result, moving toward or 
away from the rotation axis can only be achieved by sig- 
nificantly increasing the angular momentum of ascending 
fluid parcels or decreasing the angular momentum of de- 
scending fluid parcels. The timescale for this angular 
momentum exchange is longer than the typical convec- 
tion timescale for plumes to traverse the domain, and 




Fig. 2. — A pedagogical illustration, using two anelastic mod- 
els, of the importance of rotation in the brown-dwarf parameter 
regime. The left column shows a rapidly rotating model (10 hour 
rotation period), and the right column shows a slowly rotating 
model (2000 hour rotation period). A constant heat flux is ap- 
plied at the bottom boundary, leading to convection. Both models 
are fully three-dimensional simulations extending 360° in longitude 
and adopt Jovian-like radial profiles of density and thermal expan- 
sivity from the SCVH equation of state. In each model, the top 
panel shows the transient initial stage soon after convection initi- 
ates, and the bottom panel shows the state after the convection 
is well developed. Colorscale denotes entropy perturbations at an 
arbitrary longitude, shown in the radius-latitude plane. Rotation 
vector points upward in the figure. In the slowly rotating case, 
rotation plays no role in the dynamics, whereas in the rapidly ro- 
tating case, the rotation forces the large-scale flow to align along 
columns parallel to the rotation axis. 

thus convection toward or away from the rotation axis is 
inhibited. Therefore, as expected from Section [2j rota- 
tion imposes on the flow a columnar structure. 

Rotation strongly affects the vertical convective veloc- 
ities as well. To order of magnitude, convective velocities 
w and temperature perturbations 6T relate to the con- 
vective heat flux F as 



pWCp 



ST 



(21) 



where Cp is specific heat at constant pressure and p 
is the local density. Convective temperature perturba- 
tions relate to convective density perturbations Sp via 
a ST Sp/p, where a is thermal expansivity. The stan- 
dard nonrotating mixing-length scaling results from as- 
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Fig. 3. — Vertical (i.e., radial) velocities for brown dwarf con- 
vection models showing that rotation significantly affects the con- 
vective velocities. Each symbol shows the domain- averaged, mass- 
weighted root-mean-square vertical wind speed versus the mass- 
weighted mean buoyancy flux agF/pCp (units m'^s"'') for a given 
numerical integration. Circles show slowly rotating models (rota- 
tion period 100 hours), while triangles show rapidly rotating mod- 
els (rotation period 10 hours). For each rotation period, models 
with a range of buoyancy flu xes were perf ormed. The dotted and 
dashed lines show Equations \22\ and (|23|, respectively. 



suming that buoyancy forces g5p/p cause free accelera- 
tion of convective plumes over a mixi ng length I, yielding 
(e.g., |Clayton|p68l |Stevenson||1979| ) 



(agFl\ 
V PCp ) 



1/3 



(22) 



In contrast, in a rapidly rotating convective flow, convec- 
tive buoyancy forces often approximately balance vertical 
Coriolis forces. Assuming the turbulent motions are ap- 
proximately isotropic (i.e., horizontal eddy velocities are 
comparable to vertical convective velocities), one instead 
obtains a vertical velocity 
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(agF\ 

\pcpnj 



1/2 



(23) 



where we have introduced a dimensionless prefactor 7 
that is expected to be order unity. Laboratory experi- 
ments in rotating tanks demonstrate that this expression 
works well in explainin g the convective velocities in the 
rapidly ro tating regime (Golitsyn'1980||1981||Boubnov & 
Golitsyn|[T990; Fernando et al. 1991). A siinilar expres- 



sion has also been suggested for the dynamo-generating 
region of planetary interiors where a three-way force bal- 



prevail (IStarchenko & Jones 2002 


Stevenson 2003 


2010 


Showman et al. (2011) showed that it also provides 



ditions. 

Our models demonstrate that, under ty pical brown 



dwarf conditions, the rotating scaling ( 23 ) provides a 
significantly better match to the convective velocities 
than the non-rotating scaling (|22| ). This is illustrated 
in Fig. |3j which shows vertical velocities for our fully 




Latitude 

Fig. 4. — Characteristic vertical convective velocities versus lat- 
itude in a convection model with a rotation period of 10 hours, 
indicating that convective velocities tend to be greater at mid-to- 
high latitudes than at low latitudes. Velocities are calculated at a 
given latitude and pressure as the root-mean-square in longitude 
and time. 

equilibrated brown dwarf models}^ Symbols depict the 
mass-weighted, global-mean vertical velocities for a se- 
quence of models that are slowly rotating (circles, rota- 
tion period 100 hours) and rapidly rotating (triangles, 
rotation period 10 hours). They are plotted against the 
mass-weighted, global-mean buoyancy flux, agF/ pcp, for 
each model. 

There are several points to note in Figure |3] First, 
the velocities in the rapidly rotating models are smaller 
than in the slowly rotating models, indicating the rota- 
tional suppression of convective motions. Second, the 
models show that the dependence of vertical velocity on 
buoyancy flux is weaker in the slowly rotating case than 
in the rapidly rotating case. The dotted line shows the 
nonrotating scaling (22), with a mixing length of 8 16 km, 



while the dashed line gives the rotating scaling ( 23 ), with 
a prefactor 7 = 0.75. The agreement is good, showing 
that, in our slowly rotating models, the mass-weighted 
mean vertical velocities scale approxiately as buoyancy 
flux to the one-third power, whereas in our rapidly ro- 
tating models, the mass-weighted mean vertical veloci- 
ties scale approximately as buoyancy fl ux t o the one-half 
power — ^just as predicted by Equations (22) and (23), re- 
spectively. At a given buoyancy flux, the two scalings 
shown in Figure [3] differ by only a factor of a few, but 
the discrepancy becomes greater with decreasing buoy- 
ancy flux, and the two predictions differ signiflcantly for 
conditions inside a typical L/T dwarf. Adopting pa- 
rameters appropriate for an L/T dwarf (a « 10~^ K, 
Cp « lO^Jkg-iK-i, and F ~ lO'^-lO^Wm-^ corre- 
sponding to effective temperatures of ^650 to 1150 K) 
yields buoyancy fluxes appropriate to the bulk interior 
of ~10~^ to 10~^m^s~'^. For these values, the convec- 
tive velocities predicted by the rotating scaling are an 
order of magnitude lower than those predicted by the 
non-rotating scaling. 
Under conditions appropriate to a typical T dwarf. 

These, and all subsequent, models are fully equilibrated brown 
dwarf models forced by internally distributed heating and cooling. 



thereby (unlike Figure [21 avoiding 
lower thermal boundaryTayers. 



the generation of unrealistic 
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Equation (23) predicts velocities of ^O.lnis ^ in the 



deep interior, ^10ms~ at 1000 bar and ^40 ms^ at 
100 bar. However, the equation likely overpredicts the 
velocities near the top of the convection zone. In par- 
ticular, because aj p is large near the outer boundary, 
the buoyancy forces are large, and this likely implies a 
breakdown of Equation (23) in the outermost part of 
the convection zone. Interactions of convection with the 
radiative-convective boundary may also be important in 
modifying the convective velocities there, an effect not 
included in Equation ( 23 ) . 



The convective velocities tend to be greater at high 
latitudes than at low latitudes, as expected from angular 
momentum constraints. This is illustrated in Fig. H] for 
a model with a rapid (10 hour) rotation period. Radial 
convective motion at low latitudes can only occur if fluid 
parcels gain or lose significant angular momentum as 
they change distance from the rotation axis; in contrast, 
convective motion near the poles involves comparatively 
modest changes in distance from the rotation axis and 
can occur more readily. The result is greater convective 
velocities near the poles than the equator. Nevertheless, 
rotational constraints still influence high-latitude convec- 
tion: the continuity equation demands that the vertical 
convective motion must necessarily be accompanied by 
horizontal convergence and divergence and thus motions 
toward or away from the rotation axis. For convective 
motions at sufficiently large scale that Ro < 1, such ro- 
tational constraints will still play an inhibiting role in 
the efficiency of polar convection. This may explain why 
the velocities in Fig. |4] vary by only a factor of ~2 from 
equator to pole. 

We now examine the large-scale flow in our fully equi- 
librated brown dwarf models. The temperatures develop 
significant organization, as shown at the 1-bar level in 
Fig. [5] for a typical model. Over a wide range of con- 
ditions, the temperatures exhibit significant latitudinal 
gradients, with polar temperatures exceeding equatorial 
temperatures by typically a few K. This equator-to-pole 
temperature difference results from the greater efficiency 
of convection in polar regions than equatorial regions. 
Cooling to space continually decreases the entropy of 
fluid near the ~lbar level; in polar regions, this low- 
entropy fluid readily sinks and is replaced with higher- 
entropy material rising from below. But in equatorial 
regions, the suppression of radial convection inhibits the 
dense, low-entropy fluid at ^Ibar from readily sinking. 
The characteristic hot-poles-cold-equator pattern seen in 
Fig.[5]is the result. We emphasize that this effect emerges 
naturally from the dynamics and is not the result of 
any latitudinally varying forcing; indeed, our forcing and 
boundary conditions are independent of latitude. 

The temperature patterns also develop significant vari- 
ations in both longitude and latitude on regional scales of 
typically ~10'' m (Fig. [5]). Convection produces regional- 
scale thermal anomalies that vary substantially in time. 
At high latitudes, these regional anomalies tend to ex- 
hibit comparable longitudinal and latitudinal scales, as 
might be expected from the fact that the plumes move 
nearly vertically and converge or diverge horizontally in 
a quasi-isotropic fashion there. At low latitudes, how- 
ever, the structures exhibit north-south elongation. This 
is the manifestation of the columnar structure taken by 




-2 







Fig. 5. — Temperature variations at 1 bar in a typical, rapidly 
rotating brown dwarf model. Colorscale gives the temperature per- 
turbation T' at 1 bar, that is, the deviation of temperature from 
its reference value. Rotation period is 10 hours. 



the convection at relatively large scales. 

The characteristic convective temperature perturba- 
tions and equator-to-pole temperature differences in our 
models decrease with increasing rotation period. This 
is illustrated in Figure [6j The top panel shows the lon- 
gitudinal (zonal) mean temperature versus latitude at 
the 1-bar level for a model with a rotation period of 10 
hours, illustrating the hot poles and cold equator with a 
difference of ~2 K. The bottom panel shows the equator- 
to-pole temperature difference (black circles), and root- 
mean-square (rms) temperature variations (red trian- 
gles), both at the 1-bar level, for a sequence of models 
with differing rotation periods. Both the equator-to-pole 
temperature differences and rms temperature perturba- 
tions are nearly constant from rotation periods of 3 hours 
to ^^40 hours. At rotation periods exceeding ~50 hours, 
however, the temperature perturbations decrease signif- 
icantly. This results from the fact that, at long rotation 
periods, the Rossby number becomes large and the con- 
vection is no longer rotationally inhibited. 

The temperature contrasts expected in the convecting 
region can be understood by combining Eqs. (21) and 
( 23 ) to yield a relation for the convective temperature 
perturbations in a rotationally dominated flow (cf |Show^ 
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man et al.||2010[ ) 



ST 



pcpag 



(24) 



Our models are performed for Jovian-like internal pro- 
files, corresponding to = 1.74 x 10^"* K^^, Cp = 
1.3 X 10^Jkg~^K~^, and gravity, density, and thermal 
expansivity at the 1-bar level of 23ms~^, 0.2kgm~'^, 
and 0.006 K~^, respec tively. As discussed in detail by 



Showman et al. (2011 1, global convective models of giant 
planets must, tor computational reasons, be overforced 
by several orders of magnitude; our model adopts a heat 
flux near 1-bar that is close to 10^ Wm^'^. For these 
values, Eq. (24 1 predicts 5T ^ 2K, very similar to the 
values actually occurri ng i n our models (e.g.. Fig. [s]). 
This indicates that Eq. ( 24 ) provides a reasonable repre- 
sentation of the model behavior. Extrapolating now to 
the conditions of a typical L/T dwarf, we adopt a tem- 
perature of 1000 K, corresponding to a radiated IR flux 
of ~ 6 X lO'' Wm~^. Inserting parameters appropri- 
ate to the 1-bar level of a brown dwarf (p = 0.03 kg m^^, 
Cp = 1.3x 10* Jkg-iR-i, a = lO-^K, = 3x lO^^s-i, 
and g = 200 ms"^), 

we obtain ST ^ 0.5 K as the ex- 
pected convective temperature perturbation for a typical 
brown dwarf. 

The convective structure exhibits significant tempo- 
ral variability, as can be seen in Figure [7j The figure 
shows the 1-bar temperature structure (at the top of the 
convection zone) at 4.8-hour intervals in a brown dwarf 
model with a rotation period of 10 hours. Such convec- 
tive variability should cause significant variability in the 
overlying atmosphere, helping to explain the variab ility 
in lightcurve shapes observed in sever al L/T dwarfs ( |Ar-| 
tigau et al.||2009[|Radigan et al.|2012[ ). We return to tEe 
dynamics of the stratified atmosphere in Section [4] 

3.3. Large-scale flow organization in the convection 

zone 

Here we address the question of whether the convection 
zone can develop organized, large-scale horizontal winds 
such as fast east-west (zonal) jets, since these might play 
a role in causing differential zonal motion of cloudy and 
cloud-free regions in the overlying atmosphere. 

At pressures > 1 Mbar, hydrogen metallizes and 
mag netohydrod ynamic (MHD) e ffects become important 



Weir et al. 1996: Nellis et al. 1995, 1996; NcUis 2000, 



2006 1. 'I'heoretical arguments and numerical simulations 
of convection in electrically conducting spherical shells — 
as applied to Jupiter, Earth's outer core, and related 
systems — suggest that the Lorentz force acts to brake 
the large-scale east-west (zonal) winds when the electri- 



metallic region (e.f 


r.. Kirk & Stevenson|1987 


Grote et al. 


2000 


Busse 


2002 


Liu et al. 2008p. INumericai simula- 



tions of dynamo generation m convecting, rotating fluids 
at high electrical conductivity have led to scaling laws 
for the magnetic field strength of ra pidly rot ating plan- 
ets and convective stars fc hristcns en fc Au bcrt 2006; 

2009[ [CEristensen||2010| . Application 
mg laws to brown dwarfs predict 



Christensen et al 

of these these scaling laws to brown dwarfs pre dict that 
brown dwarfs will exh ibit strong magnetic fields ( Reiners 
& Christensen 20101. These dynamo experiments also 
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Fig. 6. — Top: Longitudinal-mean temperature versus latitude at 
1 bar for a model with a rotation period of 10 hours, illustrating the 
emergence of a systematic equator-to-polo temperature difference. 
Bottom: Black dots show longitudinally averaged pole-to-equator 
temperature differences, and red triangles denote the root-mean- 
square horizontal temperature fluctuations, both at 1 bar, for a 
sequence of otherwise identical models varying in rotation period 
from 3 to 200 hours. 

lead to scaling laws for th e mean flow velocities in the 



dynamo-generatin g region (Christensen & Aubert 2006_ 
Christensen|2010 )■ When heat fiuxes, rotation rates, anc 
densities appropriate to typical L/T-transition dwarfs 
are adopted (F ~ lO^Wm-^, Vl - lO^^-io-^ g-i, and 
p ~ 1-5 X lO^'kgm"'^), these scaling laws predict typi- 
cal fluid velocites of ~0.1-0.3ms~^ — similar to estimates 



from Equation (23) under the same assumptions. Over- 
all, these results suggest that the flow speeds are weak 
in the metallic interiors of brown dwarfs. We for now 
proceed under the assumption that the large-scale hor- 
izontal winds are weak in the metallic region, and ask 
what happens in the overlying molecular envelope. 

The emergence of large-scale, organized horizontal 
temperature gradients (cf Figures [5] and [6 1 implies that 
the flow will develop large-scale shear of the zonal wind 
in the direction along the rotation axis via the thermal- 
wind equation ([?]). We here write this in the form 



2n 



ghct — ~ gkjctaST 



(25) 



dz ' p 

where Sp and ST are the characteristic large-scale hori 
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zontal density and temperature differences (on isobars) 
which occur over a horizontal wavenumber fcjot (between 
the equator and pole, for example). We envision that 
these horizontal density and temperature differences re- 
sult from large-scale organization of the convective tem- 
perature f luct uations, and we therefor e eq uate 6T in 
Equation (25 1 to that from Equation (24 1. Doing so 



yields a characteristic variation of the zonal wind along 
z of 



Au ; 



fcjct 



f fFga\ 

J \pcpnj 



1/2 



dz 



(26) 



The quantity a/p in the integrand of (26) varies by or- 
ders of magnitude from the atmosphere to the deep in- 
terior and must be accounted for. In contrast, F, g, 
and Cp vary radially by a f actor of two or le ss across 



the molecular envelope (see Kaspi et al. 2009), and to 
a first approximation — here seeking simply an order-of- 
magnitude expression — we can treat them as constant. 
If we furthermore adopt the ideal-gas equation of state, 
which is reasonably accurate in the outermost layers, and 
assume that the background thermal profile is an adia- 
bat, we can integrate Equation (26) analytically to ob- 



tain the characteristic difference m zonal wind (along z) 
between a deep pressure pbot and some low pressure p: 



Al 



FR 



1/2 



1 



1 



{1 -2k)\ sin (j}\p^ ycpgfl/ 

L ^hot -I 

. (27) 

where we have used the fact that the pressure variation 
along z is overwhelmingly dom inated by the hydrostatic 
contribution. In Equation ([27]), 6 — T{po/p)'^ is the po- 
tential temperature of the adiabat (that is, a representa- 
tion of the entropy of the adiabat), po is a reference pres- 
sure (which we take here to be 1 bar), R is the specific gas 
constant, k — R/cp, and the region under consideration 
has a characteristic latitude ip. 

Adopting values appropriate to a typical brown dwarf 
{R = 3700Jkg-iK-i, K = 2/7, 9 = lOOOK, n « 
3 X 10-*s-i, (j) « 30°, F ~ lO^Wm-^ g « SOOms^^ 
and fcjot = 1 X 10~^ m corresponding to a wavelength of 
approximately one Jupiter radius), the equation can be 
expressed 



A-u ; 



1 bar 
P 



1 bar 

Pbot 



(28) 



We are interested in the wind shear between the deep 
interior and the top of the convection zone, where the 
pressure is approximately p ~ 1 bar. Interestingly, when 
we consider a ny d eep pressure pbot ^ the second term 
in Equation ( 28 1 drops out and the equation becomes 
independent oFpbot; this is because the factor a/p be- 
comes extremely small at high pressure, so that almost 
all of the contribution to Am comes from the outermost 
few scale heights of the convection zone — even if a very 
deep layer is being considered. With the adopted param- 
eters, we then obtain Au ~ 2ms~^ for the difference in 
zonal wind (in the direction of z) between any deep level 
and 1 bar. The implication is that, given the expected 
temperature variations associated with convection, the 
large-scale wind varies by at most a few ms~^ along the 
direction of the rotation axis. If the large-scale, zonal- 
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Fig. 7. — Snapshots at different times of the temperature pertur- 
bations at 1 bar in a single model with rotation period of 10 hours. 
Temperature perturbations are deviations of temperature from the 
reference state, in K. Time separation between frames is 4.8 hours. 

mean horizontal wind is weak in the deep interior where 
MHD effects predominate, then it will also be weak near 
the top of the convection zone. If, rather than adopt- 
ing a horizontal length scale 27r/fcjot of a Jupiter ra- 
dius, we instead adopt a smaller length scale (e.g., 0.1 
Jupiter radius, appropriate to the regional-scale temper- 
ature anomalies seen in Figure pi), we then conclude that 
horizontal winds of tens of m s are possible at the top 
of the convection zone. Despite the uncertainties, these 
estimates suggest that the large-scale zonal and merid- 
ional wind speeds in the convection zone are < 10^ ms~^ 
for typical brown dwarfs. 

4. DYNAMICS OF THE STRATIFIED 
ATMOSPHERE 

We have so far emphasized the convective interior, but 
infrared spectra and light curves emerge from pressures 
that for a wide range of effective temperature lie within 
the st ratified atmosphere ov erlying the convective region 



(e.g., Burrows et al. 2006a). Understanding this strati- 



fied region is therefore crucial for understanding obser- 
vations. In particular, the horizontal temperature differ- 
ences, wind speeds, and dominant flow length scales in 
this layer will control the variability in IR lightcurves, 
and vertical mixing rates will control cloudiness and 
chemical disequilibrium. Here we outline the expected 
dynamics of this stratified layer. 
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4.1. Qualitative mechanism of atmospheric circulation 

At first glance, it is not obvious that brown dwarfs 
should exhibit significant large-scale circulations in their 
atmospheres. Because they receive no external irradi- 
ation, the temperature-pressure profiles in their strat- 
ified atmospheres are determined primarily by absorp- 
tion of upwelling IR radiation from below. Since the 
interior entropy of a brown dwarf varies little with lat- 
itude, one might therefore expect that the radiative- 
equilibrium temperature profile of the stratified atmo- 
sphere should vary little with latitude, and that — at least 
at large scales — the stratified regions will be relatively 
quiescent. This constrasts significantly from the tropo- 
spheres of most solar system planets — and hot Jupiters — 
where differential stellar heating between equator and 
pole (or day and night) leads to a thermally driven at- 
mospheric circulation. 

However, the interaction of convective turbulence with 
the stable layer on brown dwarfs will perturb the strat- 
ified layer and generate a wide spec trum of atmospheric 
waves, including gravity w aves (e.g., Goldreich & Kumar 



1990 Freytag et al. 2010[ ) and Ross by waves. In solar - 
systern atmospheres, including that of Earth, Jupiter, 
Saturn, Uranus, and Neptune, such waves generated in 
the troposphere by convection and various instabilities 
propagate upward into the stratosphere. The interac- 
tion of these waves with the mean flow — in particu- 
lar, the generation, absorption, breaking, and dissipation 
of these waves — induces a large-scale circulation in the 
stratosphere. Indeed, despite the existence of equator- 
to-pole radiative (thermal) forcing in irradiated atmo- 
spheres, this mechanical, wave-induced forcing is per- 
haps the dominant driver of the stratospheric circula- 
tion on the Earth and the giant planets (for reviews, see, 
e.g., [Andrews et al.|p87l |Shepherd||2000] [20031 |Haynes| 



2005 r In a similar way, we envision that the break 



mg, absorption, and dissipation of convectively generated 
waves will drive a large-scale circulation in the stratified 
atmospheres of brown dwarfs. 
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Fig. 8. — Schematic illustration of a wave-driven atmospheric cir- 
culation, as occurs in the stratospheres of Earth, Jupiter, Saturn, 
Uranus and Neptune and as we propose occurs at regional-to-global 
scales in the stably stratified atmospheres of brown dwarfs. Grav- 
ity and Rossby waves (orange) propagate from the convective zone 
into the atmosphere, where they brealt or dissipate and induce an 
acceleration of the mean wind, here illustrated as a vector coming 
out of the page (black symbol labeled "eddy accel"). This drives 
a horizontal wind here also represented as a vector coming out 
of the page (red symbol labeled "primary wind"). The result- 
ing deviation from geostrophy drives a weak secondary circulation 
(blue contours) in the plane perpendicular to the primary wind. In 
steady state, the eddy acceleration is balanced by a Coriolis accel- 
eration (0 symbol, representing a vector pointing into the page) 
associated with the secondary circulation. Vertical motion asso- 
ciated with the secondary circulation advects entropy vertically, 
leading to horizontal temperature contrasts (labelled "hot" and 
"cold"). These temperature contrasts are precisely those needed 
to maintain thermal-wind balance with the vertical shear of the 
primary wind. Scales are uncertain but are plausibly thousands 
to tens of thousands of km horizontally and several scale heights 
vertically. 



A variety of nonlinear interactions and feedbacks en- 
hance the ability of such wave/mean- flow interactions to 
drive an atmospheric circulation. For example, vertically 
propagating waves are preferentially absorbed near crit- 
ical layers where the background flow speed matches the 
wave speed; such absorption causes an acceleration of 
the mean flow that is spatially coherent. In Earth's at- 
mosphere, this effect allows convectively generated waves 
propagating upward from the troposphere to drive zonal 
jets in the stratosphere, a phenomenon kno wn as the 
"Qua si-Biennial Oscillation" or QBO (Baldwin et al. 
2001|). A similar phenomenon has been observed on 



Jupiter (Friedson 1999J and has been su ggested to oc- 
cur in hot stars (e.g., Rogers et al.| 2012). Likewise, the 
mixing induced by breaking Rossby waves is spatially 
inhomogeneous and naturally leads to the form ation of 
jets and vortices (e.g., Dritschel & Mclntyre'2008|) . Ideal- 
ized numerical experiments of two-dimensional and strat- 
ified, three-dimensional, rapidly rotating flows demon- 
strate that random turbulent forcing can generically lead 
to the generation of la r ge-scale vortices and jets (e.g. , 
Nozawa fc Yoden'' 1997"l [Hua ng fc Robinson 19981 [Mar 



cus et al. 112000 



Smith fe Vallis„2001, 



Smith 2004; IScott 



fc Polvani||2007l |Showman||20071 [Dritschel fc Mclntyre] 
2008| [Dritschel fc Scott[[201 1| ) . 



Regardless of the details of this forcing, the rapid ro- 
tation will dominate the physical structure of such a cir- 
culation at large scales (defined here as say > 10"^ km). 
The small Rossby numbers expected at large scales im- 
ply that this circulation will be geostrophically balanced, 
that is, pressure-gradient forces will approximately bal- 
ance Coriolis forces in the horizontal momentum equa- 
tion. Here, we adopt the primitive equations, which are 
the standard equations governing atmospheric flows in 
stably stratified atmospheres when the horizontal dimen- 
sio ns greatly exceed the vertical dimensions (for reviews . 



see 



Pedlosky|p87| |Vallis||2006| [Showman et ai:]|2010 ) 
Using log-pressure as a vertical coordinate, geostrophy 
reads 



where u and v are the east-west (zonal) and north-south 
(meridional) wind on isobars, $ is the gravitational po- 
tential on isobars, x and y are eastward and northward 
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distance, respectively, and the derivatives are taken on 
constant-pressure surfaces. Here, / = 2f2 sin (j) is the 
Coriohs parameter. When combined with local hydro- 
static balance — valid in the stratified atmosphere at hori- 
zontal scales greatly exceeding vertical ones — geostrophy 
implies that the vertical shears of the horizontal wind 
relate to the horizontal te mperature gradients via the 
thermal- wind equation (e.g., Pedlosky|1987 Vallis|2006| 



dz dy 



dz dx 



(30) 



where z = — Inp is the vertical coordinate. At large 
scales, then, the development of fast winds in the 
atmosphere — particularly if the large-scale winds in the 
convection zone are weak — depends on the ability of the 
atmosphere to maintain horizontal temperature gradi- 
ents. 

What is the nature of this wave-driven circulation? 
Figure [8] provides a schematic illustration of the dynam- 
ics. Generally, the acceleration of large-scale horizontal 
winds by waves or turbulence induces deviations from 
geostrophy, leading to a mismatch between the pressure- 
gradient and Coriolis forces in the direction perpendic- 
ular to the wind. This unbalanced force drives a so- 
called "secondary circulation" in the plane perpendicular 
to the main geostrophic flow; the Coriolis forces and en- 
tropy advection caused by this circulation act to restore 
geos trophy. This standard mechanism is w ell unde r stood 
[see [Haynes et aL (|1991|) for theory, and James (1994 



pp. 100-107) or HonoE'pKM pp. 313-327) tor brieFre 



views] and provides the dynamical link between the hor- 
izontal winds, temperature perturbations, and vertical 
velocities. 

4.2. Quantitative model of atmospheric circulation 

Although the amplitude of the wave driving in brown- 
dwarf atmospheres remains unknown, the above dynam- 
ical arguments allow us to determine the relationship be- 
tween the horizontal winds, temperature contrasts, and 
vertical velocities as a Junction of the wave-driving am- 
plitude. We here construct a simple analytic theory of 
this atmospheric circulation, treating the wave-driving 
amplitude as a free parameter. The model is approxi- 
mate and makes a number of simplifying assumptions in 
the spirit of exposing the dynamics in the simplest pos- 
sible context. As a result, the model is not expected to 
be accurate in quantitative detail. Rather, the goal is to 
broadly illustrate the types of physical processes govern- 
ing the atmospheric circulation on brown dwarfs, and to 
obtain ordcr-of-magnitude estimates for the horizontal 
temperature perturbations and wind speeds, quantities 
important in shaping the observables. 

In steady state, the momentum balance in the direc- 
tion parallel to the geostrophic wind reads, to order of 
mag nitudej^ 

fv^A (33) 

Suppose, for concreteness, that the dominant geostrophic flow 
consists of zonal jets, as exist on Jupiter. The significant zonal 
symmetry of such jets suggests decomposing the flow into zonal- 
mean and deviation (eddy) components, A = A+A' . By expanding 
the zonal momentum equation and zonally averaging, we obtain the 
Eulorian-mean equation for the evolution of the zonal-mean flow, 



where v is the horizontal flow perpendicular to the main 
geostrophic flow and A is the characteristic magnitude of 
the eddy-induced acceleration of the mean flow, due to 
breaking, absorption, or dissipation of gravity or Rossby 
waves. What this equation implies is that the wave in- 
teractions with the mean flow drive a so-called secondary 
circulation, v (meridional in the case of zonal jets, radi- 
ally toward or away from the vortex center in the case of 
a large vortex). See Figure |8] for the conceptual picture. 

The associated vertical velocity can be obtained from 
the continuity equation, which for the primitive equa- 
tions in log-pressure coordinates reads 

I ^0 (34) 

which we can approximate to order of magnitude as 



vl 



w 
AS' 



(35) 



where w is the characteristic vertical velocity (in units 
of scale heights per sec), Az is the vertical scale of the 
circulation (in units of scale heights) , and I is the dom- 
inant horizontal wavenumber of the circulation. Equa- 



tions ( 33 1 and ( 35 ) imply that 



vlA'z 



lAAz 



(36) 



Thus, wave interactions with the mean flow drive large- 
scale vertical motions. 

These wave-driven, large-scale vertical motions will ad- 
vect entropy vertically, leading to the existence of hor- 
izontal temperature variations on isobars. These are in 
fact exactly the temperature differences needed to main- 
tain the wave-driven geostrophic flow in thermal-wind 
balance. To quantify, consider the thermodynamic en- 
ergy equation, which can be expressed as 

— +Vh-VhT-zu — -— = — (37) 
ot R Cp 

where T is temperature, the horizontal velocity, is 

u, over time (adopting Cartesian geometry for simplicity) 



du du du d(u'v') sd(e "u'vo') 

— = -V m V fv ^ - ^ 

dt dy dz dy dz 



(31) 



where zu = dz/dt = —dlnp/dt is the vertical velocity in log- 
pressure coordinates. Thus, the absorption, breaking, or dissipa- 
tion of waves can drive a mean flow, u. Scaling analysis of this 
equation immediately shows that, on the righthand side, the flrst 
and second terms are both order Ro smaller than the third term. 
In steady state, then, the balance in a geostrophic flow is between 
the eddy-driven accelerations and the Coriolis force associated with 
a mean meridional circulation, i.e., 



d{u'v') (9(n'ro') 



dy 



dz 



(32) 



If the flow consists predominantly of large vortices rather than 
zonal jets, one can alternately adopt a cylindrical coordinate sys- 
tem centered on a vortex, where u is the azimuthal flow around 
the vortex, v is the radial velocity (toward/away from the vor- 
tex center), and the eddy-mean-flow composition denotes an az- 
imuthal mean around the vortex (overbars), and deviations there- 
from (primes). Azimuthally averaging the azim utha l mo meii tum 
equation then leads to relationships analogous to (jsij and | |32| l. In 
either case, t he e quation can be expressed, to order of magnitude, 
as Equation l|33|. 
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the horizontal gradient operator, H the scale height, N 
the Brunt- Vaisala frequency, q the specific heating rate, 
and Cp the specific heat. In a statistical steady state, we 
envision a primary balance between the radiative heat- 
ing/cooling (righthand side) and vertical advection (right 
term on left side). If isentrope slopes are sufficiently 
large, the horizontai mixing may also contribute via the 
term v • VT. Wc write this balance as 



implies 



V3 — — = Vft • VftT. 



R 



(38) 



We parameterize radiative heating/cooling as Newto- 
nian relaxation of the temperature toward the radiative- 
equilibrium state, expressed as g/cp = (Tcq — T)/Ti.ad, 
where T^a^^z) is the radiative-equilibrium temperature 
profile and Tj-ad is a specified radiative time constant. 
Since brown dwarfs receive no external irradiation, Teq 
is to zeroth order independent of latitude and longi- 
tude. To order-of-magnitude, the characteristic devia- 
tion of temperature from its local radiative equilibrium, 
T(z)— roq(5), is comparable to the characteristic horizon- 
tal temperature difference on isobars, AThoriz- We also 
parameterize the meridional eddy mixing as a diffusive 
process, with eddy diffusivity D. To order-of-magnitude, 
we thus have 



hori: 



R 



Trad 



DfATi 



hori: 



(39) 



The physical interpretation is that vertical advection (left 
side) attempts to increase the horizontal temperature 
contrasts, whereas radiation and meridional eddy mix- 
ing (right side) both attempt to decrease the horizon- 
tal temperature contrasts. Here, w and AThoriz refer to 
characteristic magnitudes and are defined positive. Im- 
portantly, the two terms on the right side have the same 
sign, since they both act in the same direction, namely, 
to damp temperatu re d ifferences r^l 

Substituting Eq. ( 36 ) into Eq. T39|) yields 



fR 



AT, 



horiz 




(40) 



which can readily be solved to yield an expression for the 
meridional temperature difference in terms of "known" 
parameters: 



ATI 



lAASH^N^ 



horiz 



(41) 



We have yet to use the meridional momentum balance 
(or radial momentum balance in the case of a vortex), 
and doing so will allow us to solve for the zonal wind sev- 
eral scale heights above the radiative-convective bound- 
ary. To order-of-magnitude, the thermal-wind equation 

Breaking gravity waves will cause a vertical mixing that might 
be represented as a vertical diffusi on of entropy, leading to an ad- 
ditional source term in Equation | |39| l. Only horizontal variations 
in the amplitude of this mixing will act to alter ATjuj^iz. A priori, 
it is not clear how such variations will correlate with the overturn- 
ing circul ation nor how to parameterize them in the context of 
Equation l |39| . Since our goal is to describe the dynamics of the 
wave-driven circulation in the simplest possible context, we there- 
fore forgo any inclusion of this vertical mixing term here, with the 
understanding that more realistic models of the large-scale circu- 
lation will probably have to account for it. 



fi^AThorizAg 
/ 



(42) 



where Au is the characteristic difference between the 
wind speed at the radiative-convective boundary and 
some level of interest, say at the mean IR phot osph ere. 
If the former is small as suggested in Section |3.3[ Au 
would approximately give the actual wind speed at lev- 
els a bov e the rad iativ e-convective boundary. Inserting 
Eq. (liTl) into Eq. (|42|), we obtain 



Am 



PAAi^H^N^ 



DP 



(43) 



Together, Equations (36), (41 1, and (43) provide the 
expressions we seek for the vertical velocities, horizon- 
tal temperature differences, and horizontal wind speeds 
as a function of ^, I, and parameters that are either 
known or can be estimated. The Coriolis parameter, 
/, follows directly from the rotation period. For a 
brown dwarf of a given effective temperature and grav- 
ity, ID radiative-transfer models allow est imates of the 
vertical teniperature-pressure profile (e.g., Marley et al. 



1996j,2002[ 2010 



_^ Burrows et al.||1997[ |2006E ) and hence 

777V" Since the infrared photosphere is typically 1-3 scale 



heights above the radiative-convective boundary (Bur- 
rows et al.||2006b|), values of Az ~ 1-3 are most appro- 



priate 

What sets the dominant horizontal length scale of the 
flow, represented in the above theory by the wavenumber 
n One possibility is the Rhines scale, given by (Am//3)^/^ 
where /3 is the derivative of the Coriolis parameter with 
northward distance y. This is generally the scale at which 
Rossby waves impose anisotropy on the flow, and in many 
systems, it is the energy-containing scale. The Rhines 
scale controls the j et widths on Jupiter, Saturn, Uranus, 
and Neptune (e.g., Cho & Polvani 1996) as well as in a 
wide range of numerical simulati ons of stratified, rotat- 
ing tu rbulence (for a review, see |Vasavada fc Showman 
2005 1 . For typical brown dwarf rotation rates and the 
wind speeds estimated in Figure |9j this yields I ^ 3- 
6 X lO"*" m~^, corresponding to horizontal wavelengths of 
10,000 to 20,000 km. Another possibility is that I results 
from an interaction of turbulent energy transfers with 
the strong radiative and/or frictional damping. Stratified 
flows forced at small scales tend to exhibit upscale energy 
cascades, and a competition between the rate of upscale 
energy transfer and the radiative damping timescale then 
determines the dominant length scale. This possibility is 
at present difficult to quantify. Given the uncertainties, 
we simply adopt plausible values for I here, and leave a 
detailed investigation for future work. 

It is also worthwhile expressing our solutions in terms 
of the power exerted by the waves in driving the large- 
scale circulation. The characteristic power per mass ex- 
erted by the waves in driving the large-scale circulation is 
approximately AAu, and the power per unit horizontal 
area is ApAu/g, where p is the pressure at the radiative- 
convective boundary. Defining a dimensionless efficiency 
7], corresponding to the fraction of the brown-dwarf heat 
flux that is used to drive the atmospheric circulation, we 
then have rjF ~ ApAu/g, where F is the heat flux ra- 
diated by the brown dwarf. Using this constraint, the 
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solutions become 



horiz 



Au ■ 



V P 



rigF 



1/2 



HN 



R 



1 

Trad 



1/2 



DP 
lAzHN 



1/2 ■ 



and 



f vgF 



1/2 



DP 



DV 



1/2 ■ 



1/2 



(44) 



(45) 



(46) 



\ P J HN ^Trad 

If we adopt a blackbody flux, F = aT"^, where T is the 
typical photospheric temperature and a is the Stefan- 
Boltzma nn constant, and express th e radiative time con- 
stant as ( [Showman &: Guillot||2002[ ) 

pcp 



^rad 



(47) 



then Equations ( 44 1 , ( 45 1 , and M6| can be expressed as 



functions of temperature and dimensionless wave-driving 
efficiency (here for simplicity neglecting the horizontal 
diffusion term), 



ATy 



horiz 



IT' 



Am ' 



and 



f yTcp 
V 4 

2^1/2 grp7/2^ 

pHN 



1/2 



lAzHN 



1/2 



/ ' 
2(^rc,)V^ 



(48) 



(49) 



(50) 



Noting that HN is the approximate horizontal phase 
speed o f lon g -ver tical-wavelength gravity waves, inspec- 
tion of (48)-(50l makes clear that, to within factors of 
order unity, 

• AThoriz/r is ?7^^^ times the ratio of the gravity 
wave speed to the sound speed, 

• Au over the sound speed is 77^/^ times the ratio of 
the Rossby deformation radius, Lo = AzHN/f, 
to the dominant horizontal length scale of the flow, 
and 

• -07 is r]^^'^T~^\ times the ratio of the sound speed to 
the gravity wave speed (in other words, the time for 
the flow to advect vertically over a scale height is 
Ty^i/^Trad timcs the ratio of the gravity wave speed 
to the sound speed). 

For an isothermal, ideal-gas atmosphere, HN — 
Ry/T/cp, and the ratio of HN to the sound speed is 
k(1 — k), which is ~0.4 for an H2 atmosphere with 
K = 2/7. Thus, the ratio of the gravity wave speed to 
the sound speed is order unity. These arguments imply 
that, for small wave-driving efficiencies {rj I), the frac- 
tional horizontal temperature differences will be small, 
the horizontal wind speeds will be much less than the 
sound speed, and the time for air to advect vertically 
over a scale height will be much longer than the radia- 
tive time constant. 



4.3. Application of the theory to giant planets and 
brown dwarfs 

Detailed numerical simulations of convection imping- 
ing on a stable layer will be necessary to quantify the 
value of rj, but several previous studies provide con- 
straints. In a theoretical investigation of convection in- 
teracting with an overlying isothermal radiative zone, 
Goldreich & Kumar ( 1990 ) found that the fraction of 
the convective heat flux converted into gravity waves is 
approximately the Mach number associated with the con- 
vection, which may be ~0.01 for typical brown dwarfs. 
This presumably provides an upper limit on rj since only 
a fraction of the energy converted to waves is actually 
used to drive a large-scale circulation. Rough estimates 
for the Earth's stratospherej^ indicate that waves drive 
a circulation with an efficiency rj ~ lO"'^. Information is 
limited for Jupiter but likewise suggests 77 ~ 10~^ (foot- 
note 13). While future work is clearly needed, these esti- 
mates suggest that values of r/ ranging from lO""*^ to 10~^ 
may be appropriate to brown dwarfs. 

We first test the theory on Jupiter's stratospheric cir- 
culation. Jupiter's mean stratospheric temperature pro- 
file rises from the '--^llO-K tropopause minimum at 150 
mbar to ^170 K at 1 mbar pressure. Voyager and Cassini 
observations show that, throughout this pressure range, 
the temperature and zonal wind vary on characteris- 
tic horizontal (meridional) scales of ~10'*km. On these 
scales, temperatures vary by ~3-5K at most latitudes. 



reaching 10 K at a few latitudes and pressures (Simon 



Miller et al. 20061. Analysis of these observations in- 
dicates that, froni 1-100 mbar, zonal winds are ^20- 
30ms~^ over most of the planet but reach '-^130ms~^ 
in specific latitude s trips including the equato r and a 



narrow jet at 23°N (Simon-Miller et al. 2006). Verti 



cal velocities are less certain but have been estimated 
at ~10~^ms~"'^ throughout much of the stratosphere. 



reaching speeds of ^3 



( [Moreno fc Sedano||1997L 

lb apply Equations ([48|, 
we adopt / = 2.4 x 10 s~ 
tude), R = 3700Jkg-iK-i, 



' m s at high latitudes 



and (50 1 to Jupiter, 



"(appropriate to 45° lati- 
p = 0.1 bar, Cp — 1.3 x 
10"* Jkg~i K~i, and evaluate NH using an isothermal 
background temperature profile with a temperature of 
165 K. Using a length scale of 10"' km (implying I — 
6 X 10~^m~i) and an efficiency 77 ^ 10~^ (see footnote 
13), our theory predicts AThoriz ^ 3K, Au ^ 50ms~^, 
and w ^ 1 X 10"^ s~^, which for a scale height of 20 km 
implies a vertical velocity of ^H ^ 2 x 10~^ ms~^. The 
predicted meridional temperature contrasts and horizon- 

In Earth's mid and high latitudes, upwardly propagating 
waves lead to typical accelerations of the zo nal-mean zonal wind 
of order A ~ 10~^ ms~^ in the stratosphere ( [Andrews et al. |1987[ 
|Vallis|p006| Chapter 13). Typical zonal wind speeds m these re- 
gions are Au ~ 20ms~^. Adopting a pressure p ^ 0.1 bar repre- 
sentative of the lower stratosphere, this implies a power per area 
driving stratospheric motions of JtpAu/g ^ 0.2Wm~^. Since 
Earth's global-mean radiated flux is F = 240 Wm~^, the implied 
efficiency is 77 ~ A-pAu/{Fg) ~ 10^'^. For Jupiter, observational 
diagnosis of stratospheric heating patterns imply typical eddy ac- 
celerations of the zonal-mean zonal wind exceeding ~ 10~^ms~^ 
in the lower stratosphere (West et a!r|[T992[ [Moreno fc Sedanoj 
[I997t . Given a typical stratospheric wind speed of ~20ms , 
lower stratospheric pressure of 0.1 bar, and a radiated flux of 
F ^ 14Wm~^, this again implies an efflciency rj ~ lO""^. 
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tal wind speeds match the observations reasonably well. 
The predicted vertical velocity lies close to the upper 
end of the observationally inferred range, suggesting that 
our theoretical estimate may be several times larger than 
the actual global-mean vertical velocity (a mismatch that 
may result from the crudity of our parameterization of 
radiative heating/cooling). Over and above the spe- 
cific numerical comparisons, it is worth emphasizing that 
observational analysis of Jupiter's stratospheric circula- 
tions supports the overall dynamical framework adopted 
here; the primary complication is that Jupiter (unlike 
brown dwarfs) exhibits a latitudinal gradient of stellar 
irradiation, and thus its stratosphere exhibits aspects of 
both thermal and mechanically driven circulations. De- 
spite this complication, the comparison is encouraging, 
and gives us confidence in applying the theory to brown 
dwarfs and directly imaged EGPs. 

We next apply the the or y to brown dwarfs. Figure |9] 
shows the solutions (48), (49), and (50) as a function 
of temperature and 77 for parameter values representa- 
tive of brown dwarfs. The dominant horizontal length 
scale of the flow is poorly constrained, and we here 
adopt a wavelength of 10* km, similar to the jet widths 
on Jupiter and Saturn. For simplicity, we parameter- 
ize the mean atmospheric temperature profile in the ra- 
diative zone as vertically isothermal, which implies that 
NH = Ry/Tjcp. The magnitude of the horizontal diffu- 
sivity D is likewise poorly constrained; nevertheless, the 
baroclinic eddy fluxes that D parameterizes should be 
small if isentropes are nearly flat, which will be the case 
in strongly stratified regions if horizontal temperature 
perturbations are small compared to the mean tempera- 
ture. We here set D — Q, yielding an upper limit on the 
amplitude of the horizontal temperature and winds. 

Key points are as follows. First, at the ~1000K tem- 
peratures of typical L/T dwarfs, the predicted horizon- 
tal temperature perturbations on isobars are ^^5-50 K for 
plausible efficiencies. Notably, in the — limit, the 
temperature perturbations at a given rj are independent 
of lengthscale and rotation rate. To see why, note that 
the steady state results from a balance between forcing 
(energy input) and damping (energy removal). Energy 
input is represented by 77; in the model assumed here, en- 
ergy damping results solely from loss of potential energy 
by radiative relaxation, and for a given value of T,.ad, this 
energy damping depends on AThoriz but is independent 
of the rotation rate and flow lengthscales. The temper- 
ature perturbations in Figure [9] are much less than the 
mean temperature, which, given plausible stratifications, 
suggests that isentropes are relatively flat. This provides 
tentative a posteriori justiflcation for neglecting D, as 
horizontal energy fluxes due to baroclinic eddies — which 
is what D represents — tend to be modest under such con- 
ditions. 

Second, Figure |9] suggests that large-scale, organized 
horizontal wind speeds of tens to hundreds of ms~^ 
are plausible for L/T dwarfs. The range of speeds pre- 
dicted here are similar to th ose suggested from analysis 



Temperature differenceB (K) 



of brown-dwarf variability (Artigau et al. 2009 Radi 



gan et al. 2012), which range from tens to hundreds 
of ms~^. Unlike temperature differences, the predicted 
wind speeds depend on length scale and rotation rate. 
This is most clearly seen from the thermal-wind rela- 




TemperaLure (K) 

IIorizoriLal wind speed (rri/sec) 




TemperaLure (K) 
Vertical wind speed (scale height/sec) 




TemperaLure (K) 

Fig. 9. — Characteristic horizontal temperature differences (top), 
horizontal wind speeds (middle), and vertical velocities (bot- 
tom) from our solutions for large-scale, wave-driven circulation in 
the stratified atmosphere. Plotted as a function of temperature 
and dimensionless efficiency by which the waves drive the atmo- 
spheric circulation. Plotted values adopt a horizontal wavenumber 
i = 6 X 10~^m~^ (corresponding to a horizontal wavelength of 
10* km) and an isothermal vertical background temperature pro- 
file (for which NH = R^T/cp). Other parameter values are 
R = 3700 J kg~^ , Az = 2 (implying a circulation two scale 
heights tall), and / = 5 X 10"* s~^ (corresponding to a rotation 
period of 5 hours). 



tionship ( 42 1 : for a given AThoriz , shorter length scales 



imply larger horizontal pressure-gradient forces, requir- 
ing larger wind speed to maintain geostrophy. Likewise, 
a slower rotation rate implies weaker Coriolis forces for 
a given wind speed; to balance a given pressure-gradient 
force, faster winds would be needed at slower rotation 
rates. At large scales, these winds will comprise strati- 
fied, geostrophically balanced turbulence, possibly orga- 
nized into coherent structures such as vortices or zonal 
jets. ^ 

Third, for plausible efficiencies, our theory predicts 
large-scale vertical velocities of ~10^® to 10^^ scale 
heights per second (Figure [9]). Typical times for air 
parcels to ascend or descend over a scale height are there- 
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fore ~10^-10^ s. For a typical brown dwarf scale height of 
7km, this implies vertical velocities of ~0.01-0.1 ms~^. 
Importantly, this motion is not diffusive but is coherent 
vertically, comprising the ascending/descending branches 
of a large-scale overturning circulation. We empha- 
size that these vertical motions occur within the strati- 
fied regions above the radiative-convective boundary and 
are distinct from the convection below the radiative- 
convective boundary. 

Such a wave-driven overturning circulation differs fun- 
damentally from a thermally driven circulation such as 
a Hadley cell. Because the horizontal temperature dif- 
ferences result primarily from vertical advection of the 
background vertical entropy gradient, the model implies 
that, on isobars, the ascending regions are cold and de- 
scending regions are hot. This thermally indirect con- 
figuration results from the fact that the circulation is 
mechanically rather than thermally driven. Such ther- 
mally indirect, wave-driven circulations are generic fea- 
tures of planetary upper tropospheres and stratospheres. 
On Earth, the main overturning circulation in the strato- 
sphere, the Brewer-Dobson c irculation, is just such a 



phenomenon (for reviews, see Haynes 2005 Vallis 2006 
Chapter 13). Likewise, in Jupiter s upper troposphere, 
observations indicate the existence of ascending motion 
in the cold, cloudy "zones" ai id descending niotion in the 



warnier, les s-cloudy "belts" (Gierasch et al. 1986 West 



et al.||1992 [). Of course, in the brown-dwarf context, we 
have neglected the potentially important roles of latent 
heating and radiative feedbacks involving clouds, which 
may complicate the picture. 

In summary, our theory predicts that the stratified at- 
mospheres of brown dwarfs and giant planets will gener- 
ically contain geostrophic flows accompanied by a large- 
scale overturning circulation, both driven by the ab- 
sorption of waves propagating from near the radiative- 
convective boundary. This overturning circulation com- 
prises coherent, large-scale horizontal and vertical mo- 
tions despite the stratified thermal structure. The verti- 
cal entropy advection caused by this overturning motion 
generates large-scale horizontal temperature differences. 

5. OBSERVATIONAL IMPLICATIONS 

5.1. Clouds and lightcurve variability 

Our results suggest horizontal temperature variations 
(on isobars) in the stratified atmosphere up to ^50 K, 
and even in the absence of clouds, these variations could 
cause lightcurve variability. The blackbody flux from a 
region with a photospheric temperature T + AThoriz dif- 
fers from that in a region of photospheric temperature 
T by AF ^ 4crT^ AThoriz, implying that the fractional 
broadband flux variations emitted by different regions 
of the brown dwarf, due to thermal variations alone, 
are AF/F ^ 4 AThoriz /T. For a mean temperature 
T ~ 1000 K, the fractional temperature variations are 
-0.005-0.05, and AF/F ~ 0.02-0.2. Of course, if the at- 
mospheric features are much smaller than the radius of 
the brown dwarf, these variations will largely cancel and 
produce minimal variation in the disk-integrated flux on 
rotational timescales. But if the dominant length scales 
approach a brown-dwarf radius, these results suggest 
that thermal variations alone could produce variations 
in the flux of up to a few percent on rotation timescales. 



This might be relevant to brown dwarfs where cloud con- 
densation is not expected to occur at and above the pho- 
tosphere (e.g., on Y and late T dwarfs). 

Our models suggest that, when cloud condensation lev- 
els lie in the atmosphere, patchy c louds will form that 
could lead to significant variability. Freytag et al. (20101 
predicted cloud patchiness at very small scales due to 
fluctuating gravity waves, but such small-scale fluctua- 
tions will cancel out in a disk average and thus will not 
produce lightcurve variations. In contrast, our models 
predict large-scale, vertically coherent regions of ascent 
and descent associated with geostrophically balanced, 
turbulent structures such as large vortices, with time- 
mean vertical velocities up to — O.lms^^. Importantly, 
these vertical motions occur despite the fact that the at- 
mosphere is stably stratified. Regions of ascent promote 
cloud formation, whereas regions of descent transport 
dry air from above and inhibit cloud formation — ^just 
as occurs with the cloudy zones and less-cloudy belts 
on Jupiter. Near typical photospheric pressures of ~1 
bar, the continuum particle-settling regime is approxi- 
mately valid; if particles settle following Stokes flow, then 
the timescale to settle vertically by a scale height ex- 
ceeds '^lO^-lO^ s for particle sizes less than about 1 /im ^ 
Thus, given our predicted vertical velocities, particles 
with radii less than ~1 /Ltm can remain lofted in regions of 
large-scale ascent, with regions of large-scale descent ex- 
hibiting lower cloud abundances. Such large-scale cloud 
patchiness can lead to significant lightcurve variability as 
cloudy and relatively cloud-free regions rotate into and 
out of view. 

The large-amplitud e lig htcurve variations obse rved by 
Artigau et al.| ( |2009[ ) and [Radigan et al.| ( |2012[ ) exhibit 
not only variation on rotational timescales but also qual- 
itative changes to the shape of the rotational modulation 
over many rotation periods. This motivates a considera- 
tion of the range of dynamical timescales operating in the 
atmosphere. Given our predicted vertical velocities, and 
typical horizontal scales of T ~ 10"*^ km, the horizontal 
advection time, Tadv ~ L/U ~ 10^-10^ s. For a typical 
brown-dwarf rotation period of — 10** s, this suggests that 
qualitative, order-unity changes to the circulation (e.g., 
in the locations or shapes of vortices and detailed struc- 
ture of the turbulence) occur on timescales of typically 
~10-100 rotation periods. This is qualitatively consis- 
tent with the timescales over which major changes to the 



lightc urves are seen in observations (e.g., [Artigau et al. 
20091). 



For temperatures and molecular viscosities appropriate to 
L/T dwarfs, gas kinetic effects cause a significant deviation from 
Stokes fall speeds only at pressures less than ~0.1bar, indicating 
that Stokes settling represents a reasonable approximation near 
the IR photosphere. The Stokes settling velocity is uiscttlo = 
2Apga^/9/i, where Ap is the density difference between the parti- 
cles and air, g is gravity, a is particle radius, and fi is the dynamic 
viscosity of H2 air. Equating this fall speed to the mean vertical 
wind speed leads to the condition that particles can remain 
suspended if the particle size 



(51) 



Adopting a viscosity of 2 X 10 ^Pas appropriate to H2 at 1000 
K (Ackerman & Marlcy 200lJ, Ap = 3000 kgm-^', g = SOOms"!, 
H = 5 km, and ro = 10~''-10~^ s~^ yields the condition a < 0.5- 
2 ^m for particles to remain suspended. 
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5.2. Vertical mixing 

A variety of brown dwarfs, including Gl 229B, Gl 
570D, HD 3561B, and 2MASS J0415-0935, as well as di- 
rectly imaged planets such as those orbiting HR 8799 and 
2M1207, show evidence for dis equilibriuni abundances of 
CO, CH4, and/or NH3 (e.g., 'NoU et al."1997', 'Saumon 



etaLlgOOO. 2006. 2007; Stephens et al. 2009; Leggett et aT 



2007b|al 120 08 Hinz et al. 2010[ [Barman et aLpOllalHT 
Marley et al. 2012). Interestingly, in some of these cases 



models suggest that the quenching of CO occurs in the 
radiative zone rather than the convection zone, implying 
that vertical mixing in the stratified atmosphere is nec- 
essary to explain the observations. Generally, these ID 
models parameterize the mixing as an eddy diffusion and 
tune the diffusion coefficient to provide a match between 
synthetic and observ ational IR spectra. Recently, |Viss-| 
cher & Moses (2011) showed using updated reaction ki- 
netics and estimates of quench timescales that several of 
these studies underestimated the eddy diffusivities nec- 
essary to explain the inferred CO abundances; as a case 
study, they showed that an eddy diffusivity exceeding 
^10^ m^ s^^ is necessary to explain the CO abundance 
on Gl 229B. 



Freytag et al. (2010) demonstrate one mechanism for 
mixing in the stratified atmosphere, which is the vertical 
mixing induced by breaking gravity waves. Our mod- 
els in Section H] suggest an alternate mechanism: that 
the mixing results from the large-scale, wave-driven cir- 
culation. Our models predict transport times over a 
scale height of Tvcrt ~10^-10^s for a typical L/T dwarf 
(see Figure [9|. As described in section |4] these mo- 
tions are not diffusive but comprise vertically coher- 
ent advection extending potentially across several scale 
heights. At the temperatures appropriate for typical 
L/T-dwarf atmospheres, C0^CH4 chemical intercon- 
version tim escales at low pressure exceed these transport 



times cales ( Cooper & Showman 



201 1 1 , indicating that the wave-driven circulation can in- 



2006 Visscher & Moses 



duce chemical quenching of CO and CH4. 

The key point is that waves can cause two distinct 
mechanisms of vertical transport: direct vertical mixing 
associated with the wave breaking, and vertical transport 
by a large-scale circulation driven by wave interactions 
with the mean fiow. Which is larger depends on the 
wave spectrum, the detailed structure of the horizontal 
wind, the atmospheric stratification, and other proper- 
ties. In Earth's stratosphere and lower mesosphere, ver- 
tical mixing by the large-scale (meridional) overturning 
circul ation d ominates over that caused by wa ve break- 
ing (Holton||r983; Holton fc Schoeberl] [19881 , but the 
reverse may be true at even higher altitudes, in the up- 
per mesosphere. Detailed calculations will be necessary 
to evaluate the relative transport efficiency of the two 
mechanisms in the context of brown dwarfs. 

6. CONCLUSIONS 

We have demonstrated that brown dwarfs should ex- 
hibit vigorous atmospheric circulations at regional to 
global scales, and that this circulation will shape the 
cloudiness, lightcurve variability, and vertical mixing in- 
ferred in their atmospheres. Our main findings are as 
follows: 

• Given inferred rotation rates and plausible wind 



speeds, the circulation on cool brown dwarfs at 
scales larger than thousands of km will be rotation- 
ally dominated, exhibiting a horizontal force bal- 
ance between pressure gradient and Coriolis forces 
(geostrophic balance). The rapid rotation will con- 
trol many aspects of the dynamical behavior at 
large scales. 

• Using standard vorticity and angular-momentum 
arguments, we showed that, in the convective inte- 
rior, the rapid rotation and nearly isentropic con- 
ditions will lead to velocities and convective struc- 
ture aligned — at large scales — along columns par- 
allel to the rotation axis. We presented global, 
three-dimensional, anelastic numerical simulations 
of convection in brown-dwarf interiors that confirm 
this prediction. At large scales, convective veloc- 
ities in the deep interior are affected significantly 
by the rotation, leading to a different prediction for 
the scaling of convective velocities with heat flux 
than predicted by standard (non-rotating) mixing 
length theory. 

• Our 3D convection models show that convection 
occurs more efficiently at high latitudes than at low 
latitudes, and this can lead to a systematic equator- 
to-pole temperature difference in the upper con- 
vective zone that may reach ~1 K. The convection 
develops significant regional-scale structure likely 
to affect the overlying atmosphere in an observable 
manner. Scaling arguments suggest that, due to 
Lorentz-force braking in the deep, electrically con- 
ducting regions, and due to modest thermal-wind 
shear in the outermost (molecular) part of the con- 
vection zone, the horizontal winds in the convection 
zone at large scales will be modest (probably less 
than ~10^ms"^). 

• The interaction of convection with the overlying, 
stably stratified atmosphere will generate a wealth 
of atmospheric waves, and we argue that, just as in 
the stratospheres of planets in the solar system, the 
interaction of these waves with the mean flow will 
lead to a significant atmospheric circulation at re- 
gional to global scales. At horizontal scales exceed- 
ing thousands of km, such a wave-driven circulation 
will, to good approximation, be geostrophically 
balanced and may consist of large-scale stratified 
turbulence, vortices, and/or jets. The geostrophic 
flow will be accompanied by a large-scale overturn- 
ing ( "meridional" ) circulation that acts to maintain 
geostrophy. 

• We presented a simple analytic theory of this wave- 
driven atmospheric circulation which illuminates 
the dynamical mechanisms and enables order-of- 
magnitude estimates for the characteristic hori- 
zontal temperature differences, wind speeds, and 
vertical mixing rates as a function of the ampli- 
tude of the wave driving, the mean stratification 
of the atmosphere, and other properties. For plau- 
sible wave-driving efficiencies, this theory predicts 
that, at photospheric pressures, the circulation will 
comprise horizontal temperature variations of sev- 
eral to ^50K, horizontal wind speeds of ~10- 
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300 ms ^, and vertical velocities that advect air 
over a scale height in timescales of ~ 10^-10® s. 

• These models help to explain recent evidence for at- 
mospheric circulation on cool brown dwarfs. Our 
predicted large-scale vertical motion in the strat- 
ified atmosphere would lead to cloud formation 
and mix chemical constituents vertically, helping 
to explain the quenching of CO, CH4, and NH3 
inferred on many brown dwarfs. This overturning 
circulation is spatially organized, consisting of co- 
herent, large-scale regions of ascent and descent, 
indicating the likelihood of cloud p atchiness that 
could explain lightcurye var i ability ( Artigau et al. 



2009[ [Radigan et al.||2012[ |Buenzli~ al.||2012p . 

Wind speeds of plausibly tens to hundreds of m s~^ 
would advect cloudy and cloud-free regions hori- 
zontally, leading to significant changes in the shape 
of lightcurves over typically tens of rotation periods 
or more. 



Our models are of course highly idealized and there 
remains a pressing need for detailed, three-dimensional 
simulations of the atmospheric circulation and its inter- 
action with the convective interior. Coupling of such 
models to prescriptions of radiative transfer and cloud 
formation will be necessary not only to explain the par- 
ticular shapes and amplitudes of brown-dwarf lightcurves 
as a function of wavelength but also to understand the 
role of the atmospheric circulation in influencing the L/T 
transition. Future observations and models are likely to 
lead to major new insights over the next few years. 
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